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STRESSES IN ELASTIC PLATES 
OVER FLEXIBLE SUBGRADES 


Eric Reissner* 


ABSTRACT 


Comparative studies are made of the effect flexible surface layers have 
on load carrying capacity of the combination of surface layer and of subgrade. 
The analysis is based on the representation of the subgrade by a system of 
springs or, equivalently, by a heavy liquid. 

Simple approximate formulas are obtained for maximum tension stress in 
the surface layer and maximum deflection of surface layer for loads distrib- 
uted uniformly in circular areas. An important bound is established for the 
range of applicability of these formulas and this bound also roughly deter- 
mines the range of usefulness of the surface layer. 

It is shown that in the range of usefulness of the surface layer reasonable 
non-uniformity of load distribution, and built-in tensions of plausible magni- 
tude are of secondary importance. 

It is further shown that non-linear effects due to non-infinitesimal magni- 
tude of deflections are negligible as long as deflections of the surface layer 
are no larger than about one half the thickness of the surface layer. 


INTRODUCTION 


The present report has as its object a study of various theoretical ques- 
tions which arise in connection with the load carrying capacity of foundations 
consisting of a subgrade strengthened through the presence of a relatively 
thin surface layer with better strength characteristics than the subgrade it- 
self. The surface layer may be a mat, or a concrete layer, or the subsoil 
itself, improved down to a certain depth by means of various additives. 

Two basic questions confront the investigator at the outset of studies of 
this kind. They are 

(1) what to assume in regard to the nature of the behaviour of the subgrade, 

(2) what to assume in regard to the nature of the behaviour of the flexible 
surface layer. 

It appears that in all previous work of this kind the subgrade has been 
considered as obeying either Hooke’s Law of elasticity or it has been assumed 
that the subgrade behaves like a system of parallel tension-compression 
springs. The latter assumptic 1 is equivalent to the consideration of the sub- 
grade as a heavy liquid and is therefore often referred to as the liquid-sub- 
grade assumption. 

It may be considered certain that either of the above two assumptions is at 
best an approximation to the actual behaviour of the subgrade under load. 
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Past experience indicates however that the approximation is quite good and 
corresponds fairly well to reality for subsoils of normal consistencies and 
loads of such distribution as permit the return to more or less original con- 
ditions upon removal of the loads. It is clearly one of the functions of the 
flexible surface layer to result in a distribution of load on the subsoil permit- 
ting the existence of this condition for larger total loads than in the absence 
of the surface layer. 

Regarding the nature of the behaviour of the flexible surface layer it is 
necessary to carry out somewhat more specific considerations than in the 
case of the subgrade. We may distinguish, in reference to the action of a 
flexible surface the following types of action: 

1) bending action 

2) tension-surface action 

3) transverse shear action. 

Actions (1) and (3) are always present, although often one of them, generally 
(3), is negligible compared to the other. Action (2) may or may not be 
present. 

The manner in which the flexible surface may be represented for the pur- 
pose of analysis depends essentially on the length ratio, diameter of load sur- 
face to thickness of surface layer. For values of this ratio beginning at about 
unity it is possible to consider the surface layer as a “thin plate” in a manner 
which will be explained later on, rather than in the more detailed fashion of 
a three-dimensional continuum. 


Nature of results obtained. 


The quantities which are determined in what follows are primarily the 


following two quantities: 

1) the maximum value of the pressure on the subgrade at the interface of 
subgrade and surface layer, 

2) the maximum value of the tensile stress in the surface layer. 

Determination of these two quantities is carried out in their dependence on 
what are believed to be the essential parameters which enter into the problem 

1) elasticity of the surface layer. 

2) thickness of surface layer 

3) diameter of load area 

4) amount of total load 

5) deformation characteristics of the subgrade. 

In the course of the present analysis it is found that, subject to specific 
restrictions, which it is believed are generally met with in practice, certain 
formulas which were previously used can indeed be justified on wider grounds 
than was heretofore possible. While some of these restrictions may have 
been implicit in the earlier work, at least one of them, presently to be formu- 
lated, has previously been disregarded explicitly. 

Upon due consideration the following decision has been made in regard to 
the present analyisis. Since the object is the comparative study of surface 


layers endowed with a variety of properties and since previous analyses have 


shown fairly close analogies between the “elastic” and the “liquid” subgrade 
behaviour it was felt justified to base the present work on just one assumption 
concerning subgrade behaviour. As this one assumption it was natural to 
choose the simpler one, namely the assumption of a liquid subgrade, charac- 
terized by a foundation constant k. 

The following four special problems pertaining to the various possible 
properties of the surface layer have been solved. 
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(1) Infinite elastic plate carrying a load distributed uniformly over a 


finite circular area. The plate is treated by conventional plate theory. The 
results for this case have been given previously by Pickett and co-workers 
(reference 1) and the present calculations confirm the formulas of these 
authors. The present report adds to the results of reference 1 a requirement 
that these formulas not be used beyond a certain value of one of the basic 
parameters. For values of the parameter within the recommended range it 
is found possible to simplify the formulas of reference 1 to bring out more 
clearly the relative significance of thickness changes and elastic modulus 
changes for the surface layers. 

(2) Infinite elastic plate carrying a load distributed parabolically over a 
finite circular area. Since the actual nature of the load intensity distribution 
over the assumed circular area is not in actuality known it appears appropri- 
ate to investigate the differences in plate stress and foundation pressure in- 
tensity due to changes in load distribution, for fixed total load. The relevant 
calculations are carried out for one family of significant load distributions. 
It is found that up to a certain point the results for problem 1 and the present 
problem 2 are the same for practical purposes. Significant discrepancies ap- 
pear at about the same values of the basic parameter which were found criti- 
cal for problem 1. This fact lends additional support to the restrictions ad- 
vised here in connection with the use of the analytical formulas for the plate 
with uniform load distribution. The formulas for the present problem 2 have 
been specifically obtained for this report and it is possible that they have not 
been previously published anywhere. 


(3) Infinite elastic plate subjected to uniform tension in plane of plate and 
carrying uniformly distributed lateral load. In view of the fact that the in- 


stallation of the flexible surface may build into it a state of pretension (wit- 
ness the occurrence of tensile cracks in the chemically treated surfaces, due 
to shrinkage) it seemed advisable to investigate the influence of built-in ten- 
sion on foundation pressure and plate bending stresses. An explicit solution 
was obtained for this problem and it is believed that this solution is new. It 
is found that for the amount of tension possible in chemically treated surface 
layers the effect of this tension is minor. However, the analytical solution 
obtained should be of value beyond leading to recognition of this fact since it 
permits a numerical analysis of similar problems in a unified manner, begin- 
ning with the plate without built-in tension all the way up to the case where 
this tension is the dominating factor and the flexural resistance of the plate is 
negligible. 

(4) Non-linear effects. It is known that in most problems concerning trans- 
verse deflections of flat plates the relations between loads and deflections are 
linear as long as the deflections remain less than about one half of the plate 
thickness. Since plates on a flexible subgrade apparently had not heretofore 
been checked in regard to this property such an analysis is given as part of 
the present report. It is found that here too non-linearity becomes significant 
only for displacements of half the plate thickness or more. In many or most 
problems pertaining to flexible surfaces on flexible subgrades this means that 
the ordinary linear, small-deflection theory is adequate. The present report 
goes somewhat beyond establishing this result. It also includes a formulation 
of the non-linear problem which is particularly convenient for the actual cal- 
culation of non-linear corrections to the linear theory in the event that such 
calculations might become desirable. 

In order to facilitate the use of this report the main body of it is restricted 
to the statement of the various problems, the statement of the results and to 


690-3 


i 


some discussion of these results. Detailed derivations are incorporated in 
the Appendix. 


LIST OF SYMBOLS 


w = transverse plate deflection in inches 

k = modulus of subgrade reaction in pounds per cubic inch 
E = modulus of elasticity of surface layer in psi 

H =Poisson’s ratio of surface layer material 

h = plate thickness in inches 

r = radial distance measured from center of load surface 
a = radius of circular load surface in inches 

p = load intensity function in pounds per square inch 

D = plate stiffness factor, D = Eh’/12(1 - v2) 

V2 = Laplace operator, V2f = d2f/dr2 + r~ldf/dr 
OF = foundation pressure 
Op = extreme fiber bending stresses in plate 

= symbol for Bessel functions of first kind 
Y = symbol for Bessel functions of second kind 
H = symbo) for Hankel functions 


A= fi 
p= 


u) 


e = symbols for real and imaginary part of certain Bessel functions 


v) defined in references 3 and 6 


C = pressure coefficient, C = 0p(0)/p average 
”) = parameter in parabolic load distribution 
C =Cr+ (1 - 
Cp = parameter defined through equation (11') 
o = hydrostatic two-dimensional tensile stress 
N = Gh, stress resultant 
F =a stress function occurring in non-linear theory 
= dimensionless coordinate, = r/a 


INFINITE ELASTIC PLATE CARRYING A LOAD 
DISTRIBUTED UNIFORMLY OVER A 
FINITE CIRCULAR AREA 


The differential equation for this problem is 


D + kw =p 
For the present problem the load intensity function p is given by 


Po F @ 


Oo, 


where a is the radius of the circular load surface. 
To be determined are in particular 
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(1) the foundation pressure distribution %f given by 
op = kw, (3) 
(2) the bending stresses in the outer most fibers of the plate, in radial and 
circumferential direction. These are given by 


2 
o (4a) 
rar 


2 
1 Eh 1 dw a“w 
4 = (2 +y (4b) 
ar* 


Equations (3) and (4) hold independently of the form of the pressure intensity 
function and will be used also in subsequent sections of this report. 
The solution of equations (1) and (2) is of the following form 


we (5) 


where 


pz As 


and where J;, Hg(1) and H;(1) are Bessel and Hankel functions in the usual 


and 


notation. 4% 
For the evaluation of (5) it is to be noted that V i: e ° 
Write further 


4 
A = ¢ = Y k/ . 
Then, in the notation of reference 3 


Jo r) = = Up (er, p+ (Pr, 


1 


and 


Since all these functions, with subscripts zero and one, are tabulated in ref- 
erences 3 and 6 it is quite simple to plot the displacement profile w(r) as 
given by (5) in its dependence on r. 
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It is plausible, and may be confirmed by computation, that w is largest at 
the origin r =Q Therewith the foundation pressure Of as given by (3) is also 
largest for r = 0. We shall first discuss the value of the maximum foundation 
pressure 9p(0). From (3) and (5), 


op(0) Po [2 + (%/2) taf ani? (6) 


The result of equation (6) is also given in reference 1, in the first and second 
columns in Table 6 for pressure coefficients C. Allowing for differences in 
notation the pressure coefficient C is the ratio %(0)/po, 


wf, 


A= das VT pas Yu - ft wl 


The results expressed in equations (6) may alternately be expressed in terms 
of the so-called Kelvin functions, which are real and imaginary parts of vari- 
ous types of Bessel functions with argument proportional to yi. In terms of 
these functions the pressure coefficient C is 


Csl* ker' (6'*) 


Definitions and tables of the Kelvin functions ber, bei, ker, kei and their de- 
rivatives may be found conveniently in reference 4. 
At this point it is well to consider the actual behaviour of the pressure 


coefficient C as a function of pa = 4VK7D a. From Table 6 of reference 3 or 
pages 240-241 of reference 4 the following behaviour is obtained. 


Table 1. Values of pressure coefficient C 


pa 0 0.3 0.5 0.7 1,0 1.5 2.0 o 
c 0 0.034 0.090 0.166 0.305 0.55 -787 1 


It would seem to be the object of the flexible surface layer to reduce ap- 
preciably the maximum value of the load pressure on the subsoil which would 
be existing without the surface layer. Accordingly, for the flexible surface 
layer to be significant, pressure coefficients C quite appreciably less than 
unity would seem to be required. This, according to Table 1, reduces the 


range of physically significant values of pa = ./ p 4 to values not too much 
in excess of unity. Furthermore, returning to equation (5), it is only as long 
as prspa is not more than about unity that the deflection w inside the load 
area is substantially uniform, because only then do we have that Jo Ar) does 
not differ much from unity. The foregoing considerations suggest that tech- 
nical application of flexible surface layers be restricted to values of pa about 
equal to unity. Additional reasons for this restriction will become apparent 
in the section of this report following the present section. 

Polynomial approximation. As long asa is not larger than about unity it 
is possible to clarify the meaning of the pressure coefficient C as given by 
equation (6) by introducing into (6) for ker’ pa the first few terms of its 
series representation in the vicinity of the point pa =0. Formula 824.5 of 
reference 4 leads to the following approximate expression for C, 


C (ea)? 


690-6 


Equation (6''') should not be used for values of pa greater than unity. It is 
however useful for smaller values of pa. We may write (6''') with 


4 
= V »*)/m? the follow 


ing inetructive explicit form 


x 2 F 2 


In equation (6'''') the quantities at our disposal are modulus E and thickness 
h of the surface layer. It is seen that increasing the thickness h means re- 
ducing the pressure coefficients (as could hardly be expected otherwise) and 
likewise increasing the stiffness of the surface layer material, that is mak- 
ing E larger, reduces the value of the pressure coefficient. Of course, the 
use of the simple formula (6'''') requires that 


4 4 
ges V 12(1 - *) VE $1 (7) 


but this, as stated before, should roughly cover the useful range of a flexible 
surface layer. 


Bending stresses. Explicit values for bending stresses in accordance with 
equations a) will be given here only for the center r = 0 where their values 


would be expected to be largest. Moreover, when r = 0 the values of radial 
and circumferential bending stress are the same. Introduction of w from (5) 
into (4) gives 


AMAT PO (uy 4 ay] ce) 
where for brevity 


Equation (8) may also somewhat more conveniently be expressed in terms of 
Kelvin functions. The result is 


,(0) = Po y fa kei! @a (8') 


Restricting attention again to the cases for which pa is less than about unity 
we can use the first few terms of the series (824.6) of reference 4 for kei' 
pa and approximate for pa kei' pa in formula (8‘) as follows, 


kei! a) *10@(2/@ a) 


(0) 


(8'') 


Therewith 
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Use of the simple formula (8'') is allowed as long as the inequality (7) is 
satisfied, that is, as long as 


ka 

Inspection of (8'') reveals that as the thickness h of the surface layer de- 
creases the bending stress increases. Furthermore the bending stress in- 
creases as the modulus E of the surface layer increases. 

Referring to the earlier result on the foundation pressure stress a, (0) it 
is seen that increasing h is beneficial both for fF and op while increasing E 
is beneficial for of but detrimental for 0g. This result is, as it should be, in 
agreement with data contained in Figures 1 to 4 of reference 5. The only way 
in which the present results may at this point be considered to go beyond 
these earlier results is in the simple explicit analytical formulas (6''''), 
(7), and (8'*). 

As a final remark we add that the formulas (8) for op(0) are in agreement, 
for a suitably chosen value of Poisson’s ratio v, with the results in columns 
1 and 2 in Table 8 of reference 1. The function F3 in column 2 is related to 
the present formulas (8) as follows. 


a 


Numerical example. Assumptions: 


4 2 
10 
\/ 
10° x (4.5)? 


12 


From Table 1 (equation 6'') for foundation pressure coefficient C, 
Co 2.3 
From equation (8') for bending stress op(0), 


“8 = 
96x. 


(setting 


Therewith 


As pa satisfies the restriction (7) reasonable values of op(0) and of 9p(0) 
should follow from the approximate formulas (6'''') and (8''). From (6'''' ) 
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ae 
or 
a =9 in. E = 10° psi. 
h = 4.5 in. k = 102 pci. 
4, 82 
= 
Po 
| 


and (with y = ) from (8" ) 


4 
45x 10° 2 


= 10g(2.07) 5.47 


Clearly, the numerical data following from the approximate formulas are 
quite adequate, even though one is close to the upper end of what was stated 
to be their usable range. 


Infinite Elastic Plate Carrying a Load 
Distributed Parabolically over 
a Finite Circular Area 


Differential equation (1) again applies with the load distribution p given by 


Ps (9) 


The parameter 7 may have any real value whatsoever. The range from zero 
to two is of particular interest. When 7 = 1 the load distribution is uniform. 
When 7” = 0 it rises from no load at midpoint to maximum load at the outer 
rim of the load surface. When 7) = 2 the load intensity is greatest for r = 0 
and decreases to zero at r =a. The total load 


2mrpdr =7 is independent of the value of 7 . 
The following formulas are obtained after considerable transformations 


(0) 


Po = 1+ pe ker' pa 


(10) 
(0) 


(1 -1) [ee kei'pa - ga - 


When 7) = 1 then equation (10) reduces to (6" ) and equation (11) reduces to 
(8') as it should. We may write 


Sy (0) 
— = Cp + (1-9 Og (20°) 


[cp + (2 - og] (12') 
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In both formulas the second terms on the right describe the effect of non- 
uniformity of the load distribution. In order to be sure that the results for 
uniform load distribution are significant we should, it would appear, require 
that these second terms be negligible compared with the first terms. The 
following short table supplies information on this point. 


Table 2. Relative contributions to OF and oR 


pa 0 05 1.0 2.0 
Cr 0 .09 
Ct 0 -.002 -.17 
Cp 0 .16 35 


-.07 


Inspection of Table 2 reveals that when pa is appreciably larger than unity 
then the effect of non-uniformity of the applied load can no longer be neg- 
lected. Since the actual distribution may not be uniform it follows that cal- 
culations based on uniformity of load distribution should be viewed with cau- 
tion as soon as pa is greater than about unity. 


Infinite Elastic Plate Subjected to Uniform Tension 
in the Plane of the Plate and Carrying a 
Uniformly Distributed Lateral Load 


Let @ be the two-dimensional hydrostatic state of tension in the undeflected 
surface layer. Let N = ho be the corresponding stress resultants. The effect 
of this state of tension on the resistance against lateral loads is expressed by 
modifying the differential equation for the deflection w, as follows: 


- + kw =p (12) 


Attention will be restricted to the case of a uniformly distributed lateral 


load p, as given by equation (2). 
The solution w is found to be of the following form 


where 


A= ee, (14a) 


and 


(1) 
Ry ( ~) 
= (15a) 
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(15d) 
a 
1+ (N/2DX “) 


P= Aas (140) 


For the numerical evaluation of the effect of the built-in tensile stress 
o = N/h we restrict attention to the values of the foundation pressure at the 
origin, 0;(0), and to the bending stress in the plate at the origin, 99 (0). The 
following — are obtained for — two quantities 


op 1+ Ege + “2 J (16) 


cosp 


The quantities u,, v,, U,, V, occurring in (16) and (17) are tabulated func- 
tions defined in accordance with reference 3 as follows: 


( Ps u, ( + iv, ( 2,8) 
u) = (uy - + Av, + 


In the absence of built-in tension, that is for N = 0, we have according to 
(14b) cos2@ = 0, @ = 77/4 and the results reduce to those previously obtained 
for this case (pages 7 to 13). The solution (13) and (15) is valid in the form 


stated as long as|cos 2¢ |<}, that is as long as N<2ykD. We shall find by 
numerical examples that this relation is satisfied for all practically possible 
values of N. However, in the event that we wished to deduce from our result 
the limiting case of a “membrane” with D = 0 we would have to write our re- 
sults first in a different form from the present one. 

The use of the formulas (16) and (17) and the influence of pre-tension on 
stresses will be illustrated by a numerical example as follows. 


Numerical example. Assumptions: 


02 pci. 
05 psi. 


a =9 in. k= 
h = 4.5 in. E= 


Therewith, according to (14), 


1 
1 


4 
¢ ae 9 = .96 
x 4.57 
cos 
2/10? x 210° x 4.5? 3860 
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We restrict attention to cases for which g is of the order of hundreds psi 
at most. We have then |cos 2¢| <<1 or 26 =7/2 and ¢ = 7/4. We take that 
value of ¢ just beyond 77/4 for which the functions in the expressions for OF 
and 9p are tabulated. This means that we take 


~ = 57/18 = 50° 


and 
coe = cos 100° = -cos 80° = ~.1736 


Therewith 
o = 3860 x .1736 = 670 pal. 


We need further 
ein 50° = .766 cos 50° = .643 


u,(.96,50°) = + .356 
v,(.96,50°) = +.338 
U,(.96,50°) = -.785 
¥,(.96,50°) =  +.579 


and* 


Therewith, from (16) 


Po 4 766 643 — 
which means that the effect of 670 psi pretension is to reduce the maximum 
foundation pressure by about 13% compared to the value of this pressure in 
the absence of the pretension. 
Next, from (17), setting v = 1/4, 


fi 0° 96 [356 - 579 _ 1338 - = 5.25 
10* -766 6143 


which means that the effect of pretension in this case reduces the maximum 
bending stress by about 5%. 

While further numerical examples can readily be calculated in the same 
manner as for the foregoing case it is apparent that pretension stresses of 
some hundreds of psi have only a relatively minor influence on foundation 
pressure and bending stress in a flexible surface layer in the range of dimen- 
sions and rigidities which are under consideration. 


Analysis of Non-linear Effects 
Analysis of non-linear effects is to be based on the differential equations 


for finite deflections of elastic plates. These are, for the case of rotationally 
symmetric displacements, 


+ kw = = p (19a) 


- (19d) 


*U, and V; are on page 203 of reference 3. The quantities u, and Vv, are on 


page 203 of reference 6. 
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The quantities D, Vv, w, k, p and Eh are the same as those used in the linear 
small-deflection theory. The function F represents the stresses in the mid- 
plane of the plate which are absent for practical purposes for sufficiently 
small deformations. The above equations correspond to those considered on 
pages 52-70 of reference 1 except that here horizontal subgrade reaction is 
considered negligible. 

A natural procedure of solving (19a) and (19b) consists in the following 
iterative procedure. 

(1) Solve the first of the two equations for w as if the non-linear term in- 
volving F were absent. 

(2) With the value of w so obtained go into the second equation and obtain F 
in terms of w. 

(3) With these values of F and w compute the non-linear term in the first 
differential equation. If this term, computed in this manner, is small com- 
pared with p then the non-linear effect is negligible. 

A systematic approach by means of the introduction of suitable dimension- 
less variables can furnish important qualitative information in this problem 
without performance of any numerical calculations. Let a be radius of the 
load surface and assume that a is a “representative” length in the sense that 
significant changes are assumed to take place over distances of the order a. 
Observe further that if a is Py  aaapespmcennte distance then in the linear 
theory the segeer term D42A2w is of importance and the deflection is of order 
of magnitude pa’ 4/p, just as if it were a plate of radius a supported along the 
edge, instead continuously by the subgrade. On the basis of these observa- 
tions write 

r (20a) 


(20b) 


where now f is “of order unity”. Further 
2 


Pp 


2 
vis ae? + Tas (20a) 


Because of the disposition concerning a we have that v21 is also ‘of order 
unity.” 

Since we want to obtain F in first approximation from the second differen- 
tial equation we must non-dimensionalize F in such manner that the terms 
with F and w in this equation are of the same importance. Inspection reveals 
that this can be accomplished by setting 


4 
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| 
Here 
0 
0 > 
4 
Ppa 
where 


The system of differential equations to be considered is now of the form 


4 Po, 2 


where 


Po 


(21c) 


1<§ 


Tpere. are still two parameters in these equations, ka*/D and (Eh/D) 
(Pog / p)2. We know from the earlier work concerning small deflections that 
ka 4/D should be no more than unity. (We can now say that it should be no 
larger than this in order that the flexural terms in the differential equation 
can make their contribution.) The remaining parameter (Eh/D) (pga*/D)2 
indicates whether or not non-linearity is important. If (Eh/D) (ppa4/D)2 is of 
“order unity” then non-linearity is important. However, if 


then non-linearity is not important. This condition may be expressed in 
terms of deflections. Set 


when V is the order of the deflection according to the linear theory. Further 
Eh/D = 12(1 - v2) /h2. Accordingly the above condition may be written in the 


form 52 
12(1 - 1. 


Poe 


In other words, non-linear effects are negligible, as long as the deflection V 
is small compared with the thickness h of the plate. This result however re- 
quires for its validity that ka4/D be not large compared with unity. If ka4/D 
is large compared with unity, a case not considered important here, then a 
different procedure must be used for the analysis of the problem. 

It may finally be stated that in the event that non-linear corrections to the 
linear theory might be desired an appropriate systematic a would be 
to develop f and g in powers of the parameter (Eh/D) (pga*/D)2. 


CONCLUSIONS 


Comparative studies of flexible surface layers endowed with a variety of 
properties may be carried out most simply and with errors of tolerable mag- 
nitude by representing the subgrade by a system of tension-compression 
springs. This type of subgrade has the same behaviour as a “heavy liquid” 
subgrade characterized by 2 foundation constant k. 

As long as the ratio of diameter of load surface to thickness of surface 
layer, 2a/h, is about unity or larger it is permissible to consider the surface 
layer as a “thin plate.” For values of 2a/h between about one and two, 
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transverse shear deformation in the plate is of quantitative significance. It 
is recommended that at some future date the work of the present report be 
supplemented so as to coordinate the present results with results showing 
the correction effect of transverse shear deformation. 

The present report suggests that the practical usefulness of a flexible sur- 
face layer requires satisfaction of the limit relation 


4 
fe gus (z) 


where E is the modulus of elasticity of the surface layer material and the 
other quantities are defined earlier. 

As long as the above limit relation is satisfied one has as simple formulas 
representing maximum foundation pressure of(0) and maximum tension 
(bending) stress in the surface layer %p(0) with no more than a few percent 


(0) 
(12) 


error, 


¢,,(0) 


1+v 2 b,4 
(III) 


In the foregoing two formulas the quantity v is Poisson’s ratio of the surface 
layer material and fo} is the average loading stress, the total load on the sur- 
face layer being 7a“po. (It is noted that the maximum deflection w(0) is re- 
lated to in the W = Op/k.) 

As long as the limit relation (I) is satisfied we have the following assur- 
ances, as is shown in the body of the present report, 

(1) The loading stress on the subgrade is less than half of what it would 
be if the load were applied to the subgrade without the interposition of the 
flexible surface layer. 

(2) Non-uniformity of the loading stress distribution is of secondary 
importance. 

(3) Built-in tensions of reasonable magnitude in the flexible surface layer 
are of negligible effect. 

(4) Non-linear effects due to finite deformation are negligible as long as 
the transverse deflection due to load of the surface layer is smaller than 
about one half the thickness of the surface layer. 


ACKNOWLEDGMENT 


The present paper is based on a report prepared for the Engineer 
Research and Development Laboratories at Fort Belvoir, Virginia, under 
Contract DA-44-009 ENG-1863, the final report having been accepted March, 
1954. 

The author wishes to express his appreciation for helpful discussions on 
the subject of this paper to Mr. J.H. Reynolds, Jr., at that time Project Engi- 
neer at Fort Belvoir, and to his colleagues Drs. T. W. Lambe, G.E. Murray, 
and D.W. Taylor of M.I.T.’s Soil Mechanics group. 


690-15 


REFERENCES 


“Deflections, moments and reactive pressures for concrete pavements,” 
by G. Pickett, M.E. Raville, W.C. Jones and F.J. McCormick, Kansas State 
College Bulletin No. 65, October, 1951. 


2. “On the bending of axially symmetric plates on elastic foundations,” by 
P.M. Naghdi and J.C. Rowley, Proc. lst Midwestern Conf. Solid Mech. 1953, 
pp. 119-123. 


3. “Table of the Bessel functions Yp(z) and Y,(z) for complex arguments,” by 
the Computation Laboratory of the National Bureau of Standards, Columbia 
University Press, 1950. 


4. “Tables of integrals and other mathematical data.” by H.B. Dwight, The 
McMillan Company, 1947. 


5. “Design criteria for expediently solidified flexible soil surfaces.” 
Memorandum by J.H. Reynolds, Jr., E.R.D.L., Fort Belvoir, Virginia, 
October, 1953. 


6. “Table of the Bessel functions Jg(z) and J;(z) for complex arguments,” by 
the Mathematical Tables Project of the National Bureau of Standards, 
Columbia University Press, 1947. 


7. “Analytical studies of landing mats for airfields.” by G. Pickett, Report of 
Eng. Exp. Station, Kansas State College, Manhattan, Kansas, December 
1951. 


APPENDIX 


A. Infinite elastic plate carrying load distributed parabolically. Derivations. 
The following derivations contain the case of uniformly distributed load 
as a special case. 


Differential equation 


Po {1 + 2(1 rg 4 
DV*y*w + kw (1) 


Total load 


2 
2% prar = f [1+ 2(1 - 9 ] rar 


2 2 
= 2(1-n = Wa”, 
(2) 


is independent of 7 


Transition conditions 


raza: w, w', continuous 


where 


M, = row), -D(yw) 
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Can be written in the form 
r-a: w, w', v-w, continuous. 
General solution of differential equation 


4) 


= “hom * “part 


(5) 


(Po/k) [ 2(2 -4) (r*/a*)] 


“part = (6) 
O,a<¢pr 
where 


Selection of suitable constants of integration 

Since Yo(0) = ~© we must set Ag and A, equal to zero in the region r <a. 
Furthermore, according to references 3 and 6, Jg(Ar) and Jor) are conjugate 
complex. In accordance with reference 3 


Jol Ar) = ul pr, + ivy gr, W4) 


Jo( Ar) = gr, - gr, 
where ug and vo are tabulated as functions of pr. 
In order to have w a real function we write, for r< a, 


P = 2 
a 


(8) 


where Cy and C; are conjugate complex. 

In order to obtain the proper form for w in the region a < r it is necessary 
to consider the behaviour of the Bessel functions for large values of their 
arguments. We write, in accordance with usage, 


(1) (2) 
and observe that since 
Yo( Ar) 2 Up + iVy Yo( Ar) = Uo (11) 
we have that 
nf?) (Zr) = (Ar) (lz, 
For large values of the argument 


(1) (2) 1 


and consequently 


4 4 
= i = + 
(7) 
A= Vv BD’ 5 
= 
| 
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(-10 1) LE 
Var 


-1(1-1) ££ -1-1)-LE 


(2), 4 1 
Vi: Vir (14) 


Both these functions decay exponentially for large pr, while Hy‘) (xr) and 

Ho (2) (Ar) increase exponentially. The latter functions must therefore be 
excluded from the solution of the present problem. Altogether we have then a 
as a real solution valid ina<r, 


Pp 
we (de) + (Ar) | (15) 


where Co and Ce are conjugate complex. 


Bessel function formulas 
For evaluation of transition conditions it is necessary to utilize the follow- 
ing differentiation formulas 


(br) = - Ar) (17) 


Evaluation of transition conditions (3) 
Set 


Aa = P 
With w from (9) and (15) 
+ Ty Joya) + (2-9) 4 (a) 
aJ,(a) - 4(1 -9) (1) )+C a) 
AJ, + AT, 7) = Co + C2 


(19) 
In view of (7) and (18) we have that#i2 = -u42, This fact permits considerable 


simplification of (19). Combine the first and third of these equations as fol- 
lows by suitable addition 


Combine second and fourth equations as follows by suitable addition 


20, AJ - = 20, 


- (my) = -(2 - + 


(1) 2 
From this 
where 
A = 3, fu) - yu) 


The determinant A may be simplified by means of the identity 
Jo(2)¥,(z) - Jy (8)Yo(s) = 


(see reference 3, page XXII). Write (21c) with (10), 
A = (214) 


Therewith, 


Cy = [a - - + 2(2- (22b) 


For uniform load distribution for which " = 1 these expressions simplify to 


Cy = (23d) 
Introduction of (23) into (9) and (15) leads to equations (5) in the main body of 
the report. 


Expression for center deflections 
Introduce (22) into (9) according to which 


w(0) = (po/k) + Cj +t, ] 


Expression for foundation pressure gp(0) 


4+ 2re{c, } 
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7 
or 

(21a) 
(21b) 
(21c) 


4 
where, with VJ k/D a =x, 


2(1 Ty #2) (xe ) } (22a') 
Now, from (8), (10) and (11), 
1%/4 1 1 
(ke ) = Jo(ke ) + (Ke ) 


- VoD] + [role + 


and in view of page XXIV of reference 3 


1 
(xe 


Analogously, 


) - xeik- & 


1*/4 
(xe )= uy (x, (x, 
+ (x, 7) + Dy (x, 7) 


= (ea, + iker, « | 


and therewith 


4(1 - 1-i 


-2(1 - 9) (kerr - ikeix) 


Re {c,} = + kei,x) 


+ (-ker, K + kei, x) 2(1 


Write finally (reference 4, page 190-191) 
(1/ (ker'« - kei'K) and kei (1/ J2) (ker'k +kei'x). 


Therewith ker)ke kei,* = JZ 


ker)k - kei,K = JE 
690-20 


(26) 

| 


-(l- + 4-1) kei'« - 2(1 ~ kere 
2 


xker'x +(1 - 4) -2ker } 28) 


Introducing (28) in (25) finally results in the following expression for founda- 
tion pressure of(0), writing x = pa, 
dker' (1-9) { 1 
Po =1l ¢ er erga 
+ xei! + paker' pa | (29) 
Ra 
Equation (29) is the same as equation (10) in the main body of the report 


and it is discussed there. When " = 1 then (29) reduces to (6" ) in the main 
body of the report. 


Bending stresses 
In accordance with equations (4) of the main body of the text 


= + — - ' 


2(1 w } (50) 


Equations (30) will be evaluated for r = 0. From equation (9) and in view of 
(16) and (17), as long as r<a, 


P 

P 
[c, Ar) + %73,( Ar) - | (31b) 


lim Jo(r) = (32) 


there follows from 


r 


With A a = , ha = we have from (33) and (30) the following expressions 


and 
Since 
lim 
== 
P = 
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- a 


(34a) 
P +v,2 


Equations (34) may be written in the following form 


(0) 


itv Eb _ | 


Evaluation 
With #2 = ix2, from (22a) 


2 2M 
(ies, + iker, «| 
~ (werk + sxer x) } 


= fa (net K- ker, 


- (kei, x + ker, K) - 2(1 - 


2 K 


(36) 
Introduction of (36) into (35) gives 


(1 -») [ «kei x - 8x ker'y - 
(37) 
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| 


Equation (37) is modified to writing 


=») kee 4(1 -v) 4(1 


12(1 -v 
= 


-v) 


Therewith 


-° 3(1 pakei' ga 
Po -v) V kh +(1-p[ gekei'ga - 


Sker'@e _ 38) 
Ga ¢ a 

When 1 - 7 = 0 then (38) reduces to equation (8') of the main text. Otherwise 

it becomes equation (11) of the main text. 


Series developments 
From reference 4, pages 187- “~ and 184: 


2 z x 
kei t= (Log? - +++) - za. 


1,_x° 
- =) + eee 
2 32°-36 


2 
kerx= -ya - - 
4 
1 
- ) 
2. eee 
kei'xs (logs + ) 


2 
ker'x= ( logs 
+ 5G 
Therewith 


1+ xkor'ze +f 16 32 x" (39) 
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hk 4 l-v 
6 
2 6 
- 4 +--- 
x? 1 4 
5 


i 2. 
+ (loge - 


_* +2- 1, 
25.36 6 


x - Gkerx +(xker'x t1) - 


Introduction of (39) to (41) into (29) gives the following expansion for the 
foundation pressure 


(42) 


When (1 - 7) = 0 then (42) reduces to equation (6'') of the main text. Other- 
wise it shows that even for 1 - 7 # 0 the foundation pressure behaves as if 
1 - 7 = 0, as long as A a is small enough. This is in accordance with the ap- 
propriate statement made in the main body of the report. 

We finally note that equations (8') and (11) for 7,(0)/pg may be developed 
in an entirely analogous manner but omit the details. 


B. Infinite elastic plate subjected to uniform tension and carrying a uniform 
lateral load over a circular area. Derivations. 
The following. derivations, besides giving important practical information, 
illustrate the usefulness of the numerical tables in references 3 and 6. 


Differential equation 


Po» 


Solution 
,a 


Where wy, is the general solution of the homogeneous equation (1). 


2x7 kei'x- kerxs x* 
[ +6 
2 
22-2 
1 1 1 2 “ 
+ [ocd - - (log> - y )x 
| 
(1) 
(2) 
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Determination of wy, 
For brevity the subscript “h” will be omitted. The homogeneous equation 
(1) may be written in the form 


(97 + a2)( + = 0 (3) 
where 


3? = x/D, + 2- (4) 


Then, as in the case of absent initial tension, 


k/D (6) 


but now 
2 26 = - N/D (7) 


or 00s Ha - — (8) 
2VkD 


The general solution of (3) is composed of the general solutions of the two 
equations (2 +)2)f =0 and (V2 + }2)f = 0. Accordingly, if at the same time 
we make sure that all solutions are in real form, 


w= CyJol Ar) + Gol Ar) 


+ Ar) + AP) 
Tabulated form of Bessel functions in (9) 
According to references 3 and 6, 
+19 
Jol ere ) = er. 2) +ivo(er,p) (10a) 


ip 
ere ) Ul +1V,( er?) (10d) 


Approximate value of angled? 

According to (8) when N - 0 then @ = 77/4. When N > 0 then cos2¢< 0 and 
consequently 2¢>7/2 org >7/4. The value ¢ next to ¢ = 7/4 for which 
tables exist is @ = 5 7/18 = 50°. Sample calculations are made for this value 
of 

Asymptotic behaviour 

Knowledge of asymptotic behaviour of functions Jg and Yo is required to 

select proper combination in the range r >a. Set 


Jol Ar) Ar) = (de) (ere) (122) 


Jol Ar) Ar) = Ar) = ( ere?) (220) 


nfl) Ar) ~ Ar 1 + ising] 
0 ———e 


(128) 
-idr -ig + ising] 
1 1 


(2) 

(Ar) ~ s 

Vir far (12b) 
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= 


Since sing > 0 we have that Hy (1) (Ar) decays exponentially for large r and 
Ho (2) (Ar) increases exponentially for large r. The reverse is true for 
Hg!) (Xr) and Ho (2) Gr). Accordingly, in the range r > a the proper combi- 
nations of Bessel functions which must be used are Ho!) (xr) and Ho(2) (jr). 

When r< a the singular behaviour of the function Yo for r = 0 excludes 
them from the solution of the given problem. 


General solution of differential equation (1) 
In view of the above considerations the general solution of (1) is taken in 
the form 


rga: we (p/k) [1+ Ar)+ Thay Tr) ] (13a) 
we (po/k) xr) + Buf?) (Fr) (130) 


Equations (13) are the natural extensions of equations (A9) and (A15), to 
which they reduce, for the case 7 = 1, when N = 0 and therewith ¢ = 77/4. 


Evaluation of transition conditions 
Form of transition conditions is 


rea: w, w', ( continuous (14) 
Setting again 
ra =p (A18) 


and making use of the differentiation formulas (A16) and (A17) there are ob- 
tained the following transition equations: 


In order to simplify the system (14) we use (4) and write 


Therewith (14) may be brought into the following alternate form 


CoH”) (ua) + - - = 21 
(a) - a) - + 


From (16), by subtraction, 


(144) 
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t, [2+ -t, [2+ 21 


From (17), by subtraction, 


Equations (18) and (19) may be rewritten in simpler form. Taking conjugates 
(art (1) 
of everything and observing that H}“' (mr) = (u) 


(1) 1 
YoY 2+ 


1 
ye) = 0 
Therewith 
erew 
2+ a 
«4 
mere wa, ya) + BL” 


(22) 


Further transformations are carried out as follows, with the help of (5) 
to (8) 
1 


1 2. 
2+ N/D 3° 2-2 2- 


C 


1 
ge 2c08~2p - 2icos2fein2p 


2ein - 2icos2Pfein2p 


+ icos2 + icot2p) 
26i1n2p 


Introducing (22) and (23) into (21) gives 


= fa + | 
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(20) 
(24) 


We transform Cy further, 


mea 4 
Cy + tcot2g) (3, + 1¥,) 


we a(sing - icoef)(1 + icot2f)(J, + 1Y,) 


ea [ sing + cosfcot2p + 


1 coef - ising [ 4 | 
r ein2g (uy 1) + (vy > 


ea 


cosf(u, - vy) + sinf(v, + 


+ if{cosf(v, + Uy) sinf(u, - | (25) 
Expression for foundation pressure at center 


From equation (13a) 


Combination of (26) and (25) gives the following explicit formula 


sing 
v,( ¢4,0) + U,( ¢a,9) 


which is equation (16) of the main body of the report. 
Expression for bending stress at center 
Analogy between (13a) and equation (A9) indicates that the general formula 


(A35), with” = 1, remains applicable for the case of the plate under tension. 
Accordingly, 


Eb 1 2 
In this, in view of (24), 


3 


3 
[(cos3g - sin3Pcot2p) 


+ + | (Jy 
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A (26) 


Tle a)? -eing + icosp (J, + 1Y,) 


ein2g 


1 


prey 


Combination of (28) and (29) gives the following formula 


Equation (30) may be written in the alternate form 


sing 


4(1 - yz cosp 
(31) 


and this formula is equation (17) of the main body of this report. 


a 

: 
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